arXiv: 1503.07165vl [astro-ph.EP] 24 Mar 2015 


Mon. Not. R. Astron. Soc. OOP. fTl fTTR 2002) Printed 26 March 2015 (MN file v2.2) 


GTC OSIRIS transiting exoplanet atmospheric survey: 
detection of potassium in HAT-P-lb from narrowband 
spectrophotometry* 

P. A. Wilson^’^’^t, D. K. Sing^, N. Nikolov^, A. Lecavelier des Etangs^’^, F. Pont^, 
J. J. Fortney,^ G. E. Ballester^, M. Lopez-Morales®, J.-M. Desert^, A. Vidal- 
Madjar^’^ 

^CNRS, UMR 7095, Institut d’Astrophysique de Paris, F-75014, Paris, France 

^Sorbonne Universites, UPMC Univ Paris 06, UMR 7095, Institut d’Astrophysique de Paris, F-75014, Paris, France 
^Astrophysics Group, School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL, UK 
'^Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA 

^Lunar and Planetary Laboratory, University of Arizona, Sonett Space Sciences Building, Tucson, Arizona 85721-0063, USA 
^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, Massachusetts 02138, USA 
"^Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder CO 80309 


Accepted 2015 March 23. Received 2015 March 17; in original form 2015 February 17 


ABSTRACT 

We present the detection of potassium in the atmosphere of HAT-P-lb using optical 
transit narrowband photometry. The results are obtained using the 10.4-m Gran Tele¬ 
scopic Canarias {GTC) together with the OSIRIS instrument in tunable filter imaging 
mode. We observed four transits, two at continuum wavelengths outside the potassium 
feature, at 6792 A and 8844 A, and two probing the potassium feature in the line wing 
at 7582.0 A and the line core at 7664.9 A using a 12 A filter width (R ~ 650). The 
planet-to-star radius ratios in the continuum are found to be 0.1176 ± 0.0013 

at 6792 A and 0.1168 + 0.0022 at 8844 A, significantly lower than the two ob¬ 

servations in the potassium line: Rp\/Rt= 0.1248 + 0.0014 in the line wing at 7582.0 A 
and i?pi/i?*= 0.1268 + 0.0012 in the line core at 7664.9 A. With a weighted mean of the 
observations outside the potassium feature i?pi/i?*= 0.1174 + 0.0010, the potassium is 
detected as an increase in the radius ratio of Ai?pi/i?*= 0.0073 + 0.0017 at 7582.0 A 
and Ai?pi/i?*= 0.0094 + 0.0016 at 7664.9 A (a significance of 4.3 and 6.1 cr respec¬ 
tively). We hypothesize that the strong detection of potassium is caused by a large 
scale height, which can be explained by a high-temperature at the base of the upper 
atmosphere. A lower mean molecular mass caused by the dissociation of molecular 
hydrogen into atomic hydrogen by the EUV flux from the host star may also partly 
explain the amplitude of our detection. 

Key words: planetary systems - stars: individual: HAT-P-lb - techniques: photo¬ 
metric 


1 INTRODUCTION 

Transiting exoplanets provide astronomers with a unique 
opportunity to study the dynamics, structure and compo¬ 
sition of exoplanetary atmospheres. A valuable insight into 

* Based on observations made with the Gran Telescopic Canarias 
(GTC), installed in the Spanish Observatorio del Roque de los 
Muchachos of the Institute de Astroffsica de Canarias, in the 
island of La Palma, and part of the large ESO program 182.C- 
2018. 

t E-mail: pwilson@iap.fr 


the planet’s atmosphere is gained by measuring the stel¬ 
lar light transmitted through the exoplanet atmosphere as 
a function of wavelength. The amount of light obscured by 
the transiting exoplanet is dependent on the composition of 
its atmosphere, with absorbing or scattering species letting 
less light through at specific wavelengths, causing a change 
in the exoplanet radius as a function of wavelength. Hot 
Jupiters in particular, a subgroup of exoplanets that orbit 
close to their host star and exhibit inflated radii have proven 
invaluable in the characterisation of exoplanet atmospheres. 

Hot Jupiters have thus far shown great atmospheric 
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diversity with some planets being obscu red by high alti- 
tude hazes or c louds (e.g. HD 189733b, IPont et alJ l2013l 
and WASP-12b. ISing et al.lEoi3h whilst others have shown 
relati vely clear at mospheres (e.g. HD 2 09458b, ISing et ^ 
2008I. WASP -19b. iHuitson et alJ l2013l and HAT-P-llb, 


Fraine et al ][2l3) with clear signs of alkali metals at opti¬ 


cal wavelengths and water features at near-IR wavelengths. 
Amongst the sub-group of hot Jupiters that are not dom¬ 
inated by clouds and hazes, the most dominant sources 
of opacity in the optical are p redicted to be the alka li 
metals sodium and potassium (ISeager fc Sassel^ l2000ll . 
with their strong resonance doublets at 5890, 5896 A and 
7665, 7699 A respectively, with Rayleigh scattering by H 2 
dominatin g at shorter wavelengths. Sodium was first de¬ 


tected by Charbonneau et al.l (|2002l l in the atmosphere of 


HD 209458 b using the Space Telescope Imaging Spec¬ 
trograph (STIS) on the Hubble Space Tele scope (HST). 
These observations were later confirmed by I Snellen et al.l 
(l2008ll using the High Dispersion Spectrograph at the 
Subaru telescope, with th e sodium line profile character¬ 
ized by ISing et M (l2008l ~) using the HST. Exoplanetary 


spheres of HD 189733b (iRedfield et al.l I 2 OO 8 I: Jensen et al.l 

2 OIII: iHuitson et aLlboidl. WASP-17b ( 

Wood et al.l I 2 OIII: 

Zhou & Bavlissll2012lL XO-2b llSine et al.l 

2012l'l and HAT-P- 


ID lllNiKolov et al.llzui4ll . Jj^xopianetar y potassium nas been 
detected on the hot Jupiters XO-2b (ISing et al.ll201in and 
WASP-31b dSing et al.ll2ol5ll . 

Tunable Filters (TFs) have the unique capability of hav¬ 
ing the central wavelength and filter passband tuned to a 
specific value. TFs consist of a Fabry-Perot etalon made 
up of two parallel reflective surfaces. As a result of care¬ 
fully varying the separation between the two plates, the 
filter width and central wavelength can be accurately se¬ 
lected. TF have several advantages over low resolution spec¬ 
troscopy. They provide accurate differential photometry and 
have the unique advantage that they can be tuned to wave¬ 
lengths not contaminated by strong telluric lin es such as 
the p rominent O 2 lines near 6884 and 7621 A (ICatanzarol 
Il997l b Since no diffractio n gratings are use d, TFs can also 
be much more efficient (IColon et al.ll201^ . especially for 
observing atomic absorption features. Combining this tech¬ 
nique with the 10.4 m aperture of the GTC telescope makes 
it possible to study the atmospheres of planets orbiting stars 
such as HAT-P-1 (V = 9.87) which are much fainter than 
HD 209458 (V = 7.63) and HD 189733 (V = 7.65), which 
due to their apparent brightness, and large scale heights, are 
the two best studied cases thus far. 

In this study we present the detection of potass ium 
in the atmosphere of HAT-P-lb (iBakos et al.l 120071 '). a 
1.319 Rjup exoplanet with a mass of 0.525 Mjup, and 
thus a low average density of p = 0.345 g/cm®, on a 
4.47 day cir cular orbit around a GOV star at a d istance 
of 0.055 au (iBakos et al.l 120071 : [Johnson et al.ll2008f) . HAT- 
P-IB (BD+37 4734B) is a part of a visual binary system 
with the F8 companion star BD-l-37 4734A, located more 
than 450 light years away in the constellation of Lacerta. 
The host star appears to not be very active (see § 14.21) 
llKnutson et al.l[2010l : lNikolov et al.ll2014ll . The planet HAT- 
P-lb shows signs of a modest temperature inversion layer 
llTodorov et alj I 2 OI 0 I). HST /STIS observations of HAT-P- 


lb by Nikolov et al.l (l2014l l detected the sodium doublet, 


yet found no sign of potassium nor the pressure-broadened 
wings of the two alkali metals. Instead they found a flat op¬ 
tical spectrum, and a tentative absorption enhancement at 
wavelengths longer than ~ 0.85/rm. The HST Wid e Field 
Camera 3 (WFC3) observations by IWakeford et al.l (l2013f) 
detected a significant water absor ption feature i n the 1 .4 pm 
absorption band. The results bv iNikolov et al] (120141') indi¬ 
cate a strong optical absorber at higher altitudes. 

With its low density and large radius, HAT-P-lb is an 
interesting test case for interior and atmosphere models. 
These observations are a part of our larger spectrophotomet- 
ric survey aimed at detecting and comparing atmospheric 
features in transiting hot Jupiters (ESO programme 182.C- 
2018). In § 2 we describe the observations, and in § 3 we 
describe the analysis of the transit light curves. In § 4 we 
present a discussion of the results where we compare the 
observations to previous HST observations, and conclude in 
§ 5. 


2 OBSERVATIONS 

Observations were performed using the 10.4 m GTC tele¬ 
scope located at Observatorio del Roque de los Muchachos 
of the Institute de Astroffsica de Canarias on the island of 
La Palma. Narrowband imaging was done with the OSIRIS 
instrument using the red tuneable filter (operating range of 
6510 A - 9345 A) tuned to the minimum width of 12 A. 


Instrumental setup 

The OSIRIS instrument (ICeoa et al.ll2003l . l2000l : ICeDalll998f) 
consists of a mosaic of two Marconi CCD42-82 CCDs each 
with a 2048 x 4096 pixel detector separated by a 72 pixel 
gap between them. Each pixel has a physical size of 15 pm, 
which gives a plate scale of 0'.H25. The observations were per¬ 
formed without any binning with a readout frequency of 500 
kHz and a gain of 1.46 e~/ ADU on CCD2. This setup gives 
a readout noise of 8e”. The observations tuned to 6792 A 
were performed using the whole CCD2 array. Being one of 
the first OSIRIS observations, a sub-array mode was not of¬ 
fered. Later observations have since been performed in the 
sub-array mode as it reduces the read time yielding a higher 
cadence, and because the photometry of the hrst observa¬ 
tions was not improved by including other fainter stars in 
the HAT-P-1 field. The observations tuned to 7582.0 and 
7664.9 A were windowed to 600x700 pixels. To obtain the 
highest possible resolution the smallest possible passband of 
12 A was chosen. 


Observing log 

HAT-P-lb was observed on four separate nights. For all ob¬ 
servations, the companion star HAT-P-l-B (BD-l-37 4734A) 
was used for photometric comparison due to its close ~ 
iT'proximity allowing for windowed frames to be taken. The 
stars are of similar spectral type and brightness. The frames 
were also rotated to ensure that both the target and the 
comparison star were at the same radial distance from the 
centre of the CCD, such that both objects were observed at 
the same wavelength. Due to problems with the GTC dome 
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shutter, targets with an elevation above 72 degrees are sub¬ 
ject to vignetting. Since the 22 October 2009 or the 19 July 
2011 observations crossed this 72 degree limit, reaching the 
maximum elevation of HAT-P-lb at 80 degrees, a decrease 
in the raw target flux was expected. Despite this, no signs 
of vignetting were observed. 

22 October 2009: The tunable filter was centred on the 
continuum at 6792 A. The observations began at 20:42 UT, 
and ended at 02:25 UT. Due to variable seeing, ranging 
from about 1.4 to 0.7 arc seconds, the exposure time was 
frequently adjusted to avoid saturation from 10 to 20s. Be¬ 
ing the first observations with the OSIRIS instrument in 
our program, there were still issues present concerning the 
dark current. A reduction of the data was done both with 
and without dark frames. Including the darks increased the 
photometric noise slightly and were thus not used in the 
reduction. We conclude that reducing the data using dark 
frames only introduces noise for such short exposure times. 
Hardware upgrades have since solved the issue. 

19 November 2010: The tunable Alter was centred at 
7582.0 A and 7664.9 A, resulting in near simultaneous light 
curves at two wavelengths sampling the line wing of the 
potassium feature and the core of the KI D2-line. Tuning 
between the two wavelengths took 0.1s with the cadence 
of the wavelength alteration set by the readout time which 
was 5s. The observations started late at 22:15 UT due to 
high humidity and ended at 01:29 UT. Seeing was variable 
ranging from 0.8 to 1.3 arc seconds. The exposure times were 
adjusted to avoid saturation. An hour into the observations 
light cirrus clouds were present. 

19 July 2011: The tunable filter was centred on the con¬ 
tinuum at 8844 A. The observations began at 00:00 UT and 
ended at 04:06 UT on the 19*^ of July. Seeing was very vari¬ 
able during the sequence, with values ranging from 0.8 to 
2.2 arc seconds. After 02:40 UT, problems with the primary 
mirror caused a distortions of the PSF (Point Spread Func¬ 
tion). The problem was subsequently corrected at 03:10 UT. 

26 November 2013: These observations were a repeat 
of the 19 November 2010 observations and used the same 
setup. To decrease the overheads the exposure times were 
increased from 6 to 15 seconds. Observations started late at 
22:16 (due to high humidity) and ended at 01:29 UT (lower 
elevation limit of the GTC of 25 degrees). 

2.1 Reductions 

The image reductions were made using standard iraeQ rou¬ 
tines that included bias subtraction, flat field division and 
the removal of cosmic rays. An average of 0.5 to 1.7 cosmic 
rays were detected per image, except for the 6792 A obser¬ 
vations where an average count of 58 was detected. These 
events are likely not cosmic rays, but instead due to issues 
with the dark current. For the cosmic ray removal the IRAF 
routine cosmicrays was used with a threshold of 100 and 
a flux ratio of 4. A large number of calibration frames were 
used which included at least 50 bias frames and 100 dome 

^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


flats. Images with peak ADU counts above 45 000 were re¬ 
moved to ensure the photometry was performed in the linear 
regime of the detector. 

Aperture photometry was done using the apphot pack¬ 
age in IRAF using a Gaussian centring algorithm. In order to 
ensure the best possible photometry, a large range of aper¬ 
tures were explored varying both the aperture size as well 
as the dimensions of the sky annulus until the RMS scatter 
of the light curve residuals were minimised. The residuals 
were calculated by subtracting the data points from an opti¬ 
mal light curve fit generated using the Levenberg-Marquardt 
least-square algorithm implemented in the SciP 3 j 3 software 
distribution (optimize. leastsq) and by iteratively reject¬ 
ing any data points with standard deviation > 3 cr. The 
FWHM (Full Width at Half Maximum) values were mea¬ 
sured using the psfmeasure task in IRAF. The Barycen- 
tric Dynamical Time (BJDtdb) was obtained by comput¬ 
ing an average of the open and close times in UTC recorded 
in the ima ge headers and c onver ted using the online too0 
based upon l^stman et al.l (120101 '). The differential photom¬ 
etry was performed with the target and the nearby reference 
star at the same wavelength. The resulting light curves are 
shown in Fig. [T] and Fig. O 

3 ANALYSIS 

3.1 Transit light curve fits 

The transit light curves w ere generated u si ng th e ana¬ 
lytical transit equations of iMandel fc Agoll (|2002| ). Only 
one reference star was used as both stars have to be the 
same distance from the optical centre to be observed at 
the same wavelength. The best fitting parameters together 
with their associated uncertainties were determined by a 
Markov chain Mon te Garlo algorithm (MCMG) following 
I Wilson et al.| ||2014|) . The initial starting parameters were 
from iNikolov et al.l (l2014h . Possible time dependent corre¬ 
lations with airmass, FWHM and detector position were 
explored. The number of correlation parameters were de- 
termined using the Bayesian Information Griterion (BIG) 
(ISchwarjll978l ~) defined as: 

BIG = -f fclnn, (1) 

where k is the number of free parameters and n is the num¬ 
ber of data points. The correlation parameters that min¬ 
imised the BIG were used to de-trend the light curves. This 
was for all light curves, except the half transit light curves 
at 7582 A and 7664.9 A, where we fitted a linear correlation 
with FWHM. For the half transit light curves no correlation 
parameters where used as it only increased the BIG. 

For the observations at 6792 A we found ABIC = 79, 
for 7582 A ABIC = 54, for 7671 A ABIC = 60 and for 
8844 A ABIC = 387 when fitting a linear correlation with 
FWHM using one coefficient. Additional FWHM coefficients 
with higher order polynomials did not improve the fit sig¬ 
nificantly with the BIG value increasing slightly and were 
thus rejected. The residual scatter of HAT-P-lb transit light 

^ http://scipy.org/ 

® http://astroutils.astronomy.Ohio-state.edu/time/utc2bjd.html 
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Figure 1. The 6792 A and 8844 A observations (outside the potassium feature) with the raw light curve shown in the top panel, the 
modeled light curve in the middle panel and the best fit residuals shown in the bottom panel. The systematic trends seen in the 8844 A 
light curve were well modeled with a linear FWHM fit which incorporated the changing PSF. 
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Figure 2. Light curves probing the potassium feature. The first and third columns (yellow points) corresponding to observations at 
7582.0 A, and the second and fourth column (orange points) corresponding to the core of the potassium line at 7664.9 A. The half 
transits were obtained 19 November 2010 whilst the full transits were obtained on 26 November 2013. 
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curves as a function of FWHM is shown in Fig. |3]with the 
best linear fitting line shown in black. 

All the transits had the planet-to-star radius ratio, 
Rp/Rs, central transit time Tc, a baseline normalising fac¬ 
tor, N and a slope term s as free parameters. All the transits 
apart from the half transits had an addition al free FWHM 
param eter. The fixed parameters taken from iNikolov et ^ 
(I 2 OI 4 II included the period, P = 4.46529976+ (55) days, im¬ 
pact parameter b = 0.7501 and the quadratic limb-darkening 
coefficients, wi and U 2 , which were calculate d using the AT- 
LAS stellar atmospheric modelfl following ISin 3 feoiOl) . A 
quadratic limb darkening law of the following form was used 


^ = ( 2 ) 

where /(I) is the intensity at the centre of the stellar disk, 
p = cos{6) is the angle between the line of sight and the 
emergent intensity while ui and U 2 are the limb darken¬ 
ing coefficients. The stellar and orbital parameters were also 
kept fixed with Ra = 1.174 Rq, the eccentricity e = 0 and 
the scaled semi-maj or axis a/Ra = 9.853, the same as in 
INikolov et al.l (|2014l l. 

Due to the different path lengths the light has to travel 
through the tunable filter, the observed wavelength de¬ 
creases radially outward from the center of the field of view 
(FOV). At the wavelengths where sky emissions from the OI 
and OH emission lines exist, sky rings are seen. The posi¬ 
tion of these sky rings can be used to measure wavelength 
drifts, which may occur as the filter is being tuned between 
two different wavelengths. As no sky-rings were seen in the 
data, observing a wavelength shift becomes impossible. The 
typical wavel e ngths shifts measured in a similar study by 
IWilson et al.l (120141 ') . who used the same setup, was found 
to be ~ 6 A. Compared to the filter bandwidth of 12 A we 
conclude that any wavelength shift which may have occurred 
is unlikely to have a measurable impact on the derived ra¬ 
dius ratios. 


3.2 Noise Estimate 

To estimate the presence of systematic noise (red noise, 
(Jr ) in the da t a, we applied the time-averaging procedure 
of IPont et al.l (l2006tl . In the presence of uncorrelated noise 
(white noise, (J™), the standard deviation of binned residuals 
((7]v) are expected to follow the relation 


_ / M 

~ VN V M-I 


(3) 


where (ti is the standard deviation of all the residuals, N 
is the number of successive residuals points, M the number 
of bins. For data with red noise, the standard deviation is 
expected to increase by an amount Or such that for every 
bin containing N points. 


(T r 



(4) 


^ http://kurucz.harvard.edu/grids.html 
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Figure 3. The residual scatter of the HAT-P-lb light curves as a 
function of FWHM. The subplots on the left show a linear fit to 
the observed FWHM trend whilst the subplots on the right show 
the detrended residuals. 


The red noise contribution was estimated by choosing N to 
be equal to the number of points in the transit (nr), which 
varied in accordance with the cadence of the observations. 
However, due to the small number of out of transit points 
compared to the number of data points within the transit 
itself, the functional form of (Jr as a function of the number 
of N was calculated by averaging 10 separate fits to Eq. |4] 
using values of N e [10,20]. For N > 20, the standard 
deviation starts to fluctuate as the number of binned points 
become large (see Fig. 0. 

To take into account the effect of red noise on the radius 
ratio, we rescale the photometric uncertainties of the data 
set by a re-weighting factor, j3 expressed as. 



following the procedures of I Winn et al.l |(2003). 

We compare the estimated P value to the val¬ 
ues obtained using the wavelet method described in 
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Figure 4. The residual noise scatter of the HAT-P-lb light curves 
(solid coloured lines), compared to the Gaussian noise expectation 
(dashed line). Due to large uncertainties in determining cTr for a 
large number of binned points, the red noise term was calculated 
as an average of red noise models using a range of N from fO to 
20 . 


ICarter fc WinrJ (l2009l) and find good agreement with pre¬ 
vious calculated values, although the wavelet method tends 
to have slightly smaller beta values. We adopt a conserva¬ 
tive approach and select the /3 factor calculated using the 
time-averaging procedure explained above. Using the mea¬ 
sured amount of flux within the photometry aperture, we 
estimate that amongst the various sources of noise, ~ 49 - 
69 per cent of it is due to photon noise. Co-added photomet¬ 
ric precisions of 170 ppm were attained which is sufficiently 
accurate to detect atmospheric features. 


4 RESULTS AND DISCUSSION 
4.1 The detection of potassium 

We observed four transits, two of which were outside the 
potassium feature at 6792 A (R = A/AA = 566) and 8844 A 
(R=737) and two inside the potassium feature probing the 
line wing at 7582.0 A (R=632) and the line peak at 7664.9 A 



Figure 5. Trans mission spectrum of HAT-P-lb showing the STIS 
G750L data of llNikoIov et aI.I^2014^ as grey points and all the 
GTC d ata has hexagonal points. The model transmission spec¬ 
trum bv iFortnev et all 1120081 . l2niflh (shown as a blue line) assumes 
an isothermal hydrostatic uniform abundance with a equilibrium 
temperature of 1200 K. The model has been binned into 12 A 
bins and shifted vertically to best fit the HST and GTC datasets. 


(R=639). Both the raw and detrended light curves, together 
with their best fit models, are shown in Fig. [T] and Fig. 0 
with their best fit model parameters presented in Table [T] 
A transmission spectrum of HAT-P-lb showing all the in¬ 
dividual GTC data poin t s tog ether with the STIS/G750L 
data from iNikolov et al.l (|2014l ) is shown in Fig. [5] and the 
complete optical transmission spectrum showing a weighted 
mean of the potassium probing wavelengths at 7582.0 A and 
7664.9 A is shown in Fig.[^ Comparing the weighted means 
of the potassium probing wavelengths to the weighted mean 
of the observations centred at 6792 A and 8844 A, presented 
in Table [21 the detection of potassium is evident from an in¬ 
crease in the radius ratio of 0.0073+0.0017 at 7582.0 A (4.3 cr 
significance) and Ai7pi/J?*= 0.0094 + 0.0016 at 7664.9 A 
( 6.1 CT significance). 

The large increase in the radius ratio at the potassium 
probing wavelengths could be due to an unknown system¬ 
atic possibly related to the small number of out of transit 
data points for all the transits. The half transit observations 
conducted on the 19 of November 2010, displayed a strong 
potassium detection and it was suspected this could be a 
systematic related to a correlation between the half transit 
nature and the large increase in radius ratio observed at the 
potassium probing wavelengths. This prompted a new set of 
observations to be done in 2013 at the same wavelengths and 
using the same setup. The results show a remarkable con¬ 
sistency between the observations done almost three years 
apart with observations done at 7582 A showing a difference 
of Rpi/i?* 2 oio- 2 oi 3 = —0.0001 ± 0.0029 and those done at 
7664.9 A a difference of Rpi/R* 2010—2013 — 0.0014 + 0.0033. 
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Table 1. Light curve system parameters for HAT-P-lb. 


Parameter 

6792 A 

7582.0t A 

7582.0 A 

7664.91 A 

7664.9 A 

8844 A 

To [BJDtdb] 

2455127.515(20) 

2455520.456(57) 

2456623.387(68) 

2455520.456(76) 

2456623.386(52) 

2455761.585(93) 

Rpi/R* 

0.1176 ± 0.0013 

0.1247 + 0.0022 

0.1248 + 0.0018 

0.1280 + 0.0030 

0.1266 ± 0.0013 

0.1168 + 0.0022 

Ui 

0.3205 

0.2676 

0.2676 

0.2667 

0.2667 

0.2508 

U2 

0.2884 

0.2560 

0.2560 

0.2595 

0.2595 

0.2571 


8.4 X 10-5 

1.0 X lO-'^ 

9.1 X 10-5 

1.0 X 10-"“ 

1.0 X 10-"* 

8.2 X 10-5 

CT 7* 

9.0 X 10-5 

1.2 X 10-* 

1.8 X 10-^ 

1.8 X 10-"“ 

1.1 X 10-"“ 

3.4 X 10-"“ 

P 

1.47 

1.49 

2.25 

1.98 

1.47 

4.26 


f Half transits 


Table 2. Weighted means of the light curve radius ratios for 
HAT-P-lb. 


Wavelength 

Rpi/R* 

Base (6792 A and 8844 A) 

0.1174 + 0.0010 

Wing (7582.0 A) 

0.1248 ± 0.0014 

Core (7664.9 A) 

0.1268 ± 0.0012 


to 0.0005. Stellar variability is therefore unlikely 

to be the main cause behind the observed difference between 
results from data obtained at different epochs. We therefore 
regard the potassium feature robust and not due to system- 
atics errors which could have been introduced by variability. 


The potassium detection is further strengthened when com¬ 
pared to the two GTC observations outside the feature at 
6792 A (green hexagon in Fig. and 8844 A (red hexagon 
in Fig. O. Even though the 6792 A observations lack an 
ingress, and although the the 8844 A observations suffered 
from problems with primary mirror distortion causing a 
clear correlation with FWHM, both observations are highly 
consistent with the HST data. The 12 A bandwidth centered 
at 7582 A and 7664.9 A does not include telluric lines which, 
if present, may have enhanced the absorption. This is fur¬ 
ther confirmed by looking at the images themselves where 
no signs of telluric rings were be seen. 

From th e HS T observations at low resolution, 
iNikolov et al.l (|2014| 1 detected a sodium line, but did not 
see either broad wings or a broad potassium feature. Our 
narrow-band photometric measurements are at higher reso¬ 
lution than the HST data, with the detected 12 A potassium 
feature reaching to higher altitudes than the 30 A sodium 
line. For observations at higher altitudes, pressure broaden¬ 
ing is expected to become less important, ca using a narrow 
potass ium profile. These results echo those of iHuitson et al.l 
(l2012ti who detected a narrow sodium profile in the atmo¬ 
sphere of HD 189733b only at high resolution, indicating a 
narrow sodium feature. 

4.2 The effects of stellar variability 

Stellar variability is known to alter the flux received from 
the host star, thereby directly affecting the measured tran¬ 
sit depth. Unocculted star spots may increase the transit 
depths, whereas occulted star spots may cause the transit 
depth to be underestimated. HAT-P-1 b is not considered to 
be an active star, showing low chromospheric a ctivity in the 
Call H and K lines with log(i?HK) = —4.984 dBakos et al.l 
I2OO7I : [Knutson et al.l I2OI0I 1. HST/STVS data show no de¬ 
tectable spot crossings and long term variability monitor¬ 
ing from the ground ha s resulted in a 0.05 per cent upper 
limit on the variability dNikolov et al.l[2014h corresponding 


4.3 The effects of resolution 

The HST observations of INikolov et al.l (l2014ll showed no 
detection of potassium with a 75 A bin centred at 7682.5 A, 
between the between the D1 and D2 doublets, with a de¬ 
rived radius ratio of Rpi/Rt= 0.11824 + 0.00096. Compared 
to the GTC observations centred near the D2 potassium line 
{Rpi/Ri,= 0.1268 + 0.0012 at 7664.9 A), a radius ratio differ¬ 
ence of ARpi/Ri,= 0.0086 + 0.0015 was observed. Since the 
GTC observations near the core of the D2 line at 7664.9 A 
were done at a resolution of R = X/AX = 639 (AA ~ 12 A) 
compared to the HST observations which were done at a 
resolution of R=102 (AA ~ 75 A) in the same wavelength 
domain, we explored if the effects of resolution can explain 
these apparently discrepant results. 

We calculated the changes in the planet-to-star radius 
ratio as a function of bin width and resolution for the HST 
observations by stepwise binning a high resolution sodium 
and potassium profile model centred on the D2 sodium and 
D2 potassium line. The HST line spread function was mod¬ 
eled using a box function. As the bin width was decreased 
from 75 to 12 A a radius ratio increase of Ai?pi/J?*= 0.0020, 
ARpi/Ri,= 0.0053 and Ai7pi/7?*= 0.0086 was measured for 
a 1200 K, 3200 K and 5200 K isothermal potassium line 
model respectively (see Fig. (Tjl. When compared to the ob¬ 
served radius ratio difference of ARpi/Ri,= 0.0085 + 0.0012 
between the HST and GTC observations, we found that a 
~ 5200 K isothermal model is required to explain why the 
large potassium absorption was seen in the GTC data but 
was not seen in the HST data. This temperature is however 
too hot compared to previous temperature measurements, 
which makes it highly unlikely that resolution alone is re¬ 
sponsible for the difference in the planet-to-star radius ratio 
derived from the GTC and the HST data. One possible ex¬ 
planation is that if the potassium is probed higher up in the 
atmosphere, where the temp erature is higher, this wo uld 
lead to a larger scale height (IVidal-Madiar et al.ll201l3 lbh. 
This explanation is further discussed in § ES] 
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Figure 6. Transmission spectrum of HAT-P-lb with data from iNikolov et alJ ll2014h shown as blue (STIS/SG430L) and red 
(S TIS/G750L) points with the data from the GTC shown as hexagonal points. The model transmission spectra are those presented 
in lNikolov et al] ll2014h and assume an isothermal hydrost atic uniform abundance with a equilibrium temperature of 1200 K. The ma¬ 
genta line is an isothermal model bv iBurrows et ^ ll2010li with an ’extra absorbe r’ at altitude with an opac ity of 0.03 cm^ g“^ from 
0.4 to 1.0 ^m. The brown line is a 1200 K isothermal model (without TiO/VO) bv iFortnev et ahl ll2008l . [2010t) . Both spectra have been 
binned using a 12 A bin. 


Resolution (A/A A) 

2948 1179 491 197 79 39 20 



Resolution (A/AA) 

3832 1533 639 255 102 51 26 



..I I ... ■ I 
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Figure 7. Changes in the planet-to-star radius ratio as function of bin width and resolution for a 1200 K (blue), 3200 K (orange) and 
a 5200 K (red) isothermal model measured by centring a bin on the D2 sodium (left) and a D2 potassium line core (right). The vertical 
axis are in units of radius ratios above the white light curve radius ratio. The 12, 30 and 75 A bin widths are indicated by dashed, 
dot-dashed and dotted lines. The slight increase in radius ratio seen at 12 A in the left plot and 68 A in the plot on the right is due to 
the 5.98 A and 34 A separation between the D1 and D2 doublets of the sodium and potassium lines being included in the bin. 
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4.4 The effects of potassium abundance and mean 
molecnlar weight 

The large radius ratio observed in the core of the potas¬ 
sium line at high resolution in the GTC data and not seen 
in the HST data at lower resolution, could in part be due 
to an enhanced abundance of potassium high in the atmo- 
sph ere. Jupiter and Saturn both show a n increase in met¬ 
als llAtreva et al.l[200^ : iFlasar et al.ll2005l l. and an enhance¬ 
ment of metals in the atmospheres of hot Jupiter is a well 
kno wn consequence predicted by the core accretion m odel 
fe.g. ISafronov|[l97l : lLissau^ll993l : IPollack et aL|[l996h . An 
enhancement in the abundance of potassium high in the at¬ 
mosphere can however not explain the large absorption mea¬ 
sured in the core of the line. This is because an increase in 
the potassium abundance with altitude would be required 
with an abundance difference between the core and the base 
being many orders of magnitudes apart. This is unlikely as 
the ionisation of potassium is expected to increase with alti¬ 
tude, resulting in a decrease in abundance of neutral potas¬ 
sium at high altitudes. The condensation of potassium in 
the lower atmosphere, where the temperature is lower, could 
cause a depletion of potassium at low altitude, thus serving 
to increase the relative abundance higher up. The most likely 
condensate to form would be KCl which would happen at 
a temperature of ~ 600 K and at a pressure of bar 
llMorlev et al .11201^ but this is too far away from the con¬ 
ditions of the lower atmosphere to be a viable explanation. 

The photon-dissociation by EUV flux breaks molecu¬ 
lar h ydrogen into atomi c hydrogen at pressu res < 10~® 
bars IIGarcfa Munoj[2007l : iKoskinen et al.|[2013l l and this in¬ 
creases the scale height, causing an increase in the observed 
radius ratio. The radius ratio would at most increase by a 
factor ~ 2 if the potassium line observations probed these 
pressures. Photon-dissociation can therefore only partly ac¬ 
count for the large absorption observed. 


4.5 The effects of temperature 


The apparent slop e seen toward s shor ter wavelengths in 
the HST data of iNikolov et al.l (l2014fl could b e due to 
Rayleigh scattering ( Lecavelier Des Etangs et al.l I2008al f9l. 
As the composition of the gases that affect this slope is 
not well known, the temperature is derived assuming a 
Rayleigh slope described by a cross section which is in¬ 
versely proportional to the 4*^ power of the wavelength. 
As molecular hydrogen is the m ost abundant molecu le it 
is used to derive a temperature. INikolov et al.l ll20l4l cal¬ 
culate a temperature of 1884 + 603 K for a pure H 2 atmo- 
sphere by modeling the Rayleigh slo pe using the equation of 
iLecavelier Des Etangs et al.l (l2008ah . 

Using Spitzer/IRAG secondary eclip s e pho tometry (at 
3.6, 4.5, 5.8 and 8.0 /rm) iTodorov et al.1 (I 2 OI 0 II derived an 
average dayside te mperature of 1500 + 10 0 K. In the trans¬ 
mission spectrum, lWakeford_et_ajJ 1201^ find a ~ 1000 K 
isothermal model by I Fortney et al.ll2008l . l201(]|i best fits the 
water feature detected at 1.45 /rm, whereas Nikolov et al.l 
ll2014h find a 1200 K isothermal model best represents all 
their optical and near-IR data. 

The enhanced absorption measured at the potassium 
probing wavelengths co mpared to both the GTC data and 
the previous HST data (INikolov et akllioidh could be due 


to a larger temperature at the high altitude layers where 
the potassium is being probed. Assuming a temperature of 
1200 K at 0.1 bar, a temperature of 3200l^5gg K at the al¬ 
titude where the potassium D2 lines are observed would 
give a change in radius ratio of AJ?pi/7?*= 0.0087, which 
is consistent with what is observed. The larger temper¬ 
ature would correspond to the base of the thermosphere 
that is being probed by potassium at ~ 10~® to ~ 10~® 
bars, where a large increase in temperature is not un - 
common in solar system planets (iLindal et al.l[l9^ . Il985l l. 
This transition region between the lower atmosphere and 
the outermost, hottest layers of the upper atmosphere is 
referred to as the base of the thermosphere. High up¬ 
per atmospheric temperatures are particularly releva nt to 
hot-Jupiter planets as predicted by mo dels (e.g., lYellel 
|2004 iTian et"^ 200^ Gurcfa^MunM 2007^ and obse rved in 


HD 209458b ( Sing et al.l 2008: Vidal -Madiar et al.ll2011ah 
and HD 189733b i Huitson et al.ll2012l l. For resolved alkali 
lines, the temperature can be derived from scale height mea¬ 
surements in both the wings, which probe lower altitudes, 
and the core, which probes higher altitudes. 


5 CONCLUSIONS 

We present four GrC/OSIRIS transit observations aimed 
at probing the presence of potassium in the atmosphere 
of HAT-P-lb. Two separate transit observations, observed 
three years apart, detect the potassium feature at high con¬ 
fidence providing a detection with a 4.3 cr significance at 
7582.0 A and a 6.1 a significance at 7664.9 A relative to 
the GTC measurements outside the feature at 6792 A and 
8844 A which are both consistent with the HST data. This is 
only the second exoplanet (after XO-2b) where the presence 
of both atmospheric Na and K has been detected. The potas¬ 
sium detection in XO-2b and HAT-P-lb have both been 
from the ground. We investigate the effects resolution has 
on the measurements and argue that although the effect is 
clearly present, it is not large enough by itself to account for 
the non-detection of potassium in the HST data at a lower 
resolution. An increase in potassium abundance with alti¬ 
tude can not explain the strong potassium detection as an 
increase many orders in magnitude would be required. The 
strong presence of potassium inferred from the observations 
could however be influenced by the dissociation of molecular 
hydrogen into atomic hydrogen thereby decreasing the mean 
molecular weight and increasing the scale height. 

Finally we discuss how an increase in temperature at 
the altitude where the potassium line is being probed could 
lead to the strong detection of potassium. It is likely that 
we detect potassium higher up in the atmosphere where the 
temperatures are higher and the scale height is larger. Fu¬ 
ture observations aimed at sampling the slope of the potas¬ 
sium profile, will provide more stringent constraints on the 
derived temperature at the altitude where the potassium is 
detected. Such observations will also allow the degeneracies 
amongst the possible effects responsible for the large absorp¬ 
tion in the line core to be minimised. 
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